SU(3) dibaryons in the Einstein-Skyrme model 
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SU(3) collective coordinate quantization to the regular solution of the B = 2 axially symmetric 
Einstein-Skyrme system is performed. For the symmetry breaking term, a perturbative treatment 
as well as the exact diagonalization method called Yabu-Ando approach are used. The effect of 
the gravity on the mass spectra of the SU(3) dibaryons and the symmetry breaking term is studied 
in detail. In the strong gravity limit, the symmetry breaking term significantly reduces and exact 
SU(3) flavor symmetry is recovered. 
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I. INTRODUCTION 

The Skyrme model [l| is considered as an unified the- 
ory of hadrons by incorporating baryons as topological 
solitons of pion fields, called skyrmions. The topological 
charge is identified as the baryon number B. Performing 
collective quantization for a B ^ I skyrmion, one can 
obtain the proton and the neutron within 30% error [3|. 
Correspondingly multi-skyrmion (B > I) solutions are 
expected to represent nuclei [3, [j] . 

The Einstein-Skyrme (ES) model can be thought of as 
a model of hadrons in which baryons interact with black 
holes. The early studies [^, [g, 0] have shown that the 
Schwarzschild black hole can support spherically sym- 
metric (hedgehog) Skyrme hair which is a first counter 
example to the no-hair conjecture. Axially symmetric 
regular and black hole skyrmion solutions with B = 2 
were found in Ref. I3|. Also some multi-skyrmion solu- 
tions B > 2 have been obtained [9]. 

If a skyrmion is regarded as a nucleon or a nuclei, it 
must be quantized to assign quantum numbers like spin, 
isospin, etc.. Study of gravitational effects on the quan- 
tum spectra of skyrmions was initiated in Ref. pUJ] by 
performing collective quantization of a i? = 1 gravitating 
skyrmion. The study about a B = 2 axially symmetric 
skyrmion immediately followed to it |ll| . It was shown 
there that the qualitative change in the mass difference, 
mean charge radius and densities under the strong gravi- 
tational infiuence confirms the attractive feature of grav- 
ity while the reduction of the axial coupling and tran- 
sition moments indicates the gravitational effects as a 
stabilizer of baryons. 

In Ref. |Tl|, we investigated the gravity effects to the 
SU(2) solutions of the ES model and found out that the 
effects are seen especially in the heavier dibaryon spec- 
tra. In this paper we shall study the gravity effects to 
the SU(3) dibaryons because they have much larger mass 
compared with that of SU(2) dibaryons in terms of the 
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large fiavor symmetry breaking. The effect will be more 
apparent on such dihyperons. 

The extension of the Skyrme model into SU(3) have 
extensively been studied. Roughly speaking, solutions 
of the SU(3) are classified into two categories: one 
is embedded, another is non-embedded ones. For the 
non-embedded solutions, in Ref.^fJ] the classical S0(3) 
chiral field is simply extended into three flavor space. 
The authors constructed the solutions in the ES system 
and found particle-like and black hole solutions. More 
systematic analysis in this direction has been done in 
Refs.fy, [lj|. By using the harmonic map ansatz, the 
authors have investigated various large flavor and multi 
winding number solitons. 

On the other hand, for the embedded solutions, the 
successful approaches are based on perturbations. In 
the bound state approach [l3| static baryons are con- 
sidered as bound states of the kaon and the skyrmions. 
The chiral field is expanded in the kaon fluctuations 
around the classical solutions. This approximation works 
well for the large symmetry breaking. Collective coordi- 
nate approach is essentially a natural extension of the 
SU(2) collective quantization scheme, including symme- 
try breaking terms. Good quantum numbers for spin and 
isospin are obtained by quantizing the their rotational 
zero modes [151, lla] ■ Yabu-Ando (YA) treatment ^17| is 
a sort of the collective quantization but the collective 
Hamiltonian is exactly diagonalized by using Euler-angle 
parameterization of the SU(3) rotations. As a result, 
baryon appears to be its lowest irreducible representa- 
tion (irrep) but contain significant admixture of higher 
representations. The application to the B = 2 axially 
symmetric skyrmions |18| and to the multi-skyrmions 
(3 < -B < 8) 19] in flat space have been studied. They 
observed large mixing of higher irreps especially in ex- 
cited or higher winding number states. 

In this Letter, we construct SU(3) axially symmetric 
skyrmions in the ES system in terms of the collective 
quantization and YA approach. We estimate mass spec- 
tra of the dibaryons down to strangeness S — —6. Since 
our concern is about contributions of higher irreps to 
the lowest states and we shall show change of the mix- 
ing probabilities about varying coupling strength of the 



gravity. 



Correspondingly, the following axially symmetric 
ansatz is imposed on the metric [20| 



II. CLASSICAL GRAVITATING B = 2 
SKYRMIONS 

In this section we shall introduce basic formalism and 
classical solutions of the axially symmetric gravitating 
skyrmions. The SU(3) extended Skyrme Lagrangian cou- 
pled with gravity can be written as 

L = Lg + Ls + LsB + Lwz (1) 

where Lq is the standard Einstcin-Hilbcrt Lagrangian 



Lg = / d^ry/^ 
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and the remaining parts are about the SU(3) extension 
of the Skyrme Lagrangian, which are defined in terms 
of the chiral field U{x) G SU(3) and using the notation 
/„ — U^dnll as follows: 
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where F^ and e are basic model parameters which indi- 
cate the pion decay constant, a dimensionless parameter, 
respectively. The Lsb is comprised of all chiral and fla- 
vor symmetry breaking terms. The Lwz is usual Wess- 
Zumino term which concerns with the topological charge 
of the soliton. Ff^,m^ are the kaon decay constant and 
the mass. Nc means a number of color. 

In this Letter, we employ the solution of the chiral field 
with axially symmetric, winding number two which is 
considered as a possible candidate for the B = 2 minimal 
energy configuration [3| 

Uo{r)^exp[iF{r,9)T-nR]. (6) 

where n^ is defined as 

nR = {sin Q (r, 9) cos nip, sin <9{r, 9) sin n^p, cos Q{r, 9)) {7) 

and n e Z is the winding number. We explore the so- 
lution with n = 2. SU(3) chiral field is constructed by 
trivial embedding: 



U{r) = 



Uo{r) 
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where the metric functions / , m and / are the function 
of coordinates r and 9. This metric is symmetric with 
respect to the z-axis {9 = 0). Substituting these ansatz 
to the Lagrangian ([1]), one obtains the following static 
(classical) energy for the chiral fields 
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+ "^ [Uda^FdeQ - dgFd^ef sin' F 
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where dimensionless variable x = eF-j^r and /3^ = |^ are 
introduced. 

For the profile functions, the boundary conditions at 
the a; = 0, oo are imposed 



F{Q,9) =7r, ^(00,6*) =0, 

9,6(0, 0) = a,e(oo,0) = o. 



At 6* = and 7r/2, 



deF{x, 0) = deF{x, ■ 



0, 



e(a;,0) = 0, e(.T,|) = |. 



(11) 
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For the solutions to be regular at the origin a; = and to 
be asymptotically flat at infinity, the following boundary 
conditions must be imposed 

a,/(0, 9) = a,m(0, 9) = dj{0, 9) = 0, (15) 
/(oo,6') = ra(oo,6') ^^(oo,^) = 1. (16) 

For the configuration to be axially symmetric, the follow- 
ing boundary conditions must be imposed at = and 

7r/2 



def{x, 0) = dem{x, 0) = del{x, 0) = 0, 



(17) 



de.f{x, |) = dem{x, |) = del{x, |) = 0. (18) 



The covariant topological current is defined by 

iv{U-^V ^UU-^V pUU-^V aU). (19) 
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FIG. 1; The profile functions F, O for a = 0.080 in tfie cylindrical coordinate with dimensionless variables P := ef-j^p, Z = e/^z. 



Substituting the ansatz ([6]), ([7]) into (fT9| the zeroth com- 
ponent is estimated as 



B" = - 



TT V^.9 



■.s\T?Fs\nQ{d.^FdgQ-deFd^Q). (20) 



The baryon number of the soliton B is derived from its 
spatial integral, thus 



B= I d^r^/^B^ 
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The inner and outer boundary conditions [Fq,Qq) — 
(7r,0) and (^1,81) = (0, tt) yield B = 2. 

By taking a variation of the static energy (fTO)) with 
respect to F and 0, one obtains the equations of motion 
for the profile functions. The field equations for the met- 
ric functions /, m and I are derived from the Einstein 
equations. The explicit form of the equations is essen- 
tially same (except for contribution of the mass term) as 
reported in Ref.[ll|. 

The effective coupling constant of the Einstein-Skyrme 
system is given by 



= 47rGi^i 



(22) 



III. THE SU(3) COLLECTIVE QUANTIZATION 

We study the SU(3) extension of the axially symmet- 
ric i? = 2 skyrmions by Yabu-Ando approach together 
with naive collective coordinate quantization. SU(3) chi- 
ral field is constructed by trivial embedding 



Uir,t) = Ait)[^^^f>^\]A\t) 



(23) 



where C/g is introduced in Eq.®. A{t) is time dependent 
SU(3) rotational matrix and R{t) describes a spatial ro- 
tation of the soliton. We introduce the angular velocities 
Vta , oJi which are defined by 



a=l 
3 



(24) 
(25) 



Substituting the chiral field (|23p into the Lagrangians 
(I3|)-(l5|) and after a lengthy calculation, one finally obtain 
the effective Lagrangian of the form: 



which is the only free parameter. 

To solve the equations of motion numerically, the re- 
laxation method with the typical grid size 100 x 30 are 
performed. We observe that the solution survives at 
< a < 0.120. Including the mass term, the range 
becomes a little narrow. We show examples of our nu- 
merical results for the profile functions F, O in FiglT] and 
also for the metric functions f,l,m in Fig|21 
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FIG. 2: The metric functions /, I, m in the cylindrical coordi- 
nate with dimensionless, rescaled variables ^ — P/{l + P),( — 
Z/{l + Z),P = eUp,Z^e,Uz. 



where 7^?, Jj, /a, Is are called the moments of inertia and 
their explicit forms are 
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Is = T^ I dxd6{l-cosF) 



^x^sin^ 



+^a;^ sin 6i{(9^F)' + sin^ ^(9^,6)^} 
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^7(1 — Dgs) exhibits strength of the symmetry breaking 
and the explicit form of 7 is 



27T_F„ I o2 



(/3« - pi) \ dxdQx^ sin6'(cosF - 1), (31) 
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where /3k = ™'^'^'° . I?88 is a component of Wigner func- 
tion which is defined as 



Dab{A) = iTr(A,AtAbA). 



(32) 



From (l26l) the Hamiltonian reads 



1 



1 



,2, 



1/1 1 n^\o 3n 1 

+ -( L^ B^ + ess 33) 

where eigenvalues of diagonal operators are already in- 
serted. Here, the eigenvalue of J is spin, / is isospin, N 
is right isospin derived from iV = ipg where (poi^o) is 
the minimal irrep and L is the third component of the 
body fixed spin operator which determine parity P of the 
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FIG. 3: The coupling constant dependence of the mass difference from the classical energy for the multiplets {10}, {27}, 
{35}, {28} ((p, q) = (0, 3), (2, 2), (4, 1), (6, 0), respectively) are sho-wn in the unit of Ge'V. The quantized energy is computed in 
terms of the naive perturbation technique. 
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FIG. 4: The coupling constant dependence of the mass difference from the classical energy for the multiplets {10}, {27}, 
{35}, {28} ((p, q) — (0, 3), (2, 2), (4, 1), (6, 0), respectively) are shown in the unit of Ge'V. The quantized energy is computed in 
terms of YA. 
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FIG. 5: Mixing probability of the higher representations to 
the lowest energy state of A'^A'^ channel in {10}. 



In terms of Eg. (136)) . the expectation value of the 
Casimir invariants C2{SU{3)) in Eq. ([M)) is easily ob- 
tained 

{C2{SU{3))) = i(p2 + q^+pq + 3(p + q)) . (38) 



For the symmetry breaking term /s7(l — Dg^), the es- 
timation of the expectation value can be done by per- 
forming the integral of the three Wigner rotation ma- 
trices [22, [2^ which is evaluated by the SU(3) Clebsch- 
Gordan coefficient, or the isoscalar factor 
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Computations of the Clebsch-Gordan coefficients can be 
performed by the numerical algorithm of Ref . [2J| . 



state by the relation of P = (—1)^. The eigenvalue of 
the EsB is derived from following eigenequation 

{C2[SUi3)] + Isiil - Dss)}-^ = ess* (34) 

where C2[SU{3)] is Casimir operator of SU(3). 

We shall investigate (|34p in two folds. One is to treat 
the symmetry breaking term pcrturbatively, another is 
to diagonalize the whole via a basis of the SU(3) Wigner 
functions. 



A. Perturbative method 

If symmetry breaking effect 7 is certainly small, pertur- 
bative treatment seems to be good approximation. We 
introduce a wave function of the form 12111 



B. Diagonalization of the collective Hamiltonian 

If symmetry breaking effect is crucial, the naive per- 
turbation will substantially fail. Yabu-Ando approach 
can improves the situation. In YA, state of baryon ap- 
pears to be its lowest irrep but contain large admixture 
of higher irreps. We shall see that such mixing reduces 
in large gravity limit. 

The wave function of the Hamiltonian is expanded in 
terms of a wave function of the lowest representation 
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In terms of the basis, the eigenvalue problem in Eq. 
can be reduces to a matrix diagonalization problem. 



* 



'^nsY.NNsY^,JjM)D''l-nN3{R-') (35) 



$ 



IIaY.NNaYR,JJay^) 
IIsY.NNaYrt 



v/dM(_l)^^+^3 






(36) 



where the dimension of the (p, q) irrep is expressed by 
d(™) ^ [p+ l)[q + l)(p + q + 2)/2, the m is representa- 
tion of SU(3) group, and subscript of the Y and Yr is 
hypercharge and right hypercharge respectively. 

With the operation of the collective quantization, the 
angular velocity f^s appears linear in Eq.f^H). Therefore 
we obtain a constraint 



(37) 



which means that the symmetry C/fl(l) is redundant. 
Thus we obtain Yr = 2. 



IV. NUMERICAL RESULTS 

For the actual calculations, we fix F^ = 108 MeV, 
e — 4.84, [3-^ — 0.263. The kaon decay constant, ex- 
perimentally, is F^ sa v^-FV, but for the simplicity, we 
employ F^^ = F,^. For the kaon mass, we employ the 
experimental value, i.e., (3,^ — 0.952. 

We estimate mass spectra belonging to SU(3) multi- 
plets {10}, {27}, {35}, {28} (or in the (p, q) representation, 
(0, 3), (2, 2), (4, 1), (6, 0)^spectively). The Finkelstein- 
Rubinstein constraints [23] tells us that for {10}, {35}, 
J = 1 is chosen for the ground state, otherwise one can 
set J = [26]. The eigenvalue of L concerns with the 
third component of body fixed spin operator [3|. Sub- 
stantially it is related to the orbital angular momentum 
but no experimental identification has been done. There- 
fore in our analysis we put L = Q for all multiplets states. 
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In YA treatment, one needs to truncate the base in 
finite size. We expand the collective wave function with 
iV < 3, except for the states {S,I) = (-6,0) in {28}, 
(-4,0), (-4,1) and (-5,1/2) in {35}. In those states, 
we expand the base with A^ < 4 for obtaining sufficient 
convergence. 

In Figl3] presents the a dependence of the mass spec- 
tra within the naive perturbation scheme. Actually, we 
show mass difference between the quantized mass spectra 
and the classical energy. Also FigHlpresents the results of 
Yabu-Ando treatment. In both results the spectra as well 
as their differences within each multiplct decrease mono- 
tonically with increasing a. On the other hand, mass dif- 
ferences between different multiplet but with same quan- 
tum numbers (S*, /) increase, which have been already 
observed in the calculation of SU{2) [11|. One easily 
observe that difference between the results of two treat- 
ments disappears with increasing a. This behavior is eas- 
ily understood. In Figl5]we illustrate the mixing proba- 
bility of the multiplet for NN channel in {10} for various 
a. For increasing a, the mixing of higher representations 
are significantly decreased; as a result, the naive pertur- 
bation is sufhcient for the analysis of SU(3) dibaryons for 
such strong gravity region. In our analysis, the pion and 
the kaon mass and the coupling constants are fixed by 
their experimental values (we simply set F^, = F-^ for the 
coupling constant) and if we take into account variations 
of the mcsonic data about change of gravity, exact SU(3) 
flavor symmetry will attain at a strong gravity limit. 



stant dependence of the energy spectra of the SU(3) 
dibaryons. We have used the collective quantization in 
three flavor space. To treat the symmetry breaking term, 
we employ the lowest order (naive) perturbation to that 
as well as Yabu-Ando treatment. Both treatments have 
shown that mass differences between spectra with dif- 
ferent strangeness decrease monotonically and increase 
within different multiplet but with same quantum num- 
bers {S, I) with increasing a. In the strong gravity limit, 
the SU(3) flavor symmetry recovers; all the spectra de- 
generate in each multiplet and no mixing between the 
multiplets occur. Such symmetry restorations may be 
observed in high energy experiment at LHC. 

In this Letter, we treat a as a free parameter. In the 
Einstein-Skyrme theory, the Planck mass is related to 
the pion decay constant F^^ and coupling constant a by 
Mpi — Ft^^Ai: ja. To realize the realistic value of the 
Planck mass, the coupling constant should be extremely 
small with a ~ 0(10"'^^). However, we have shown that 
the effects of gravity can be observed only in large a. 
Some theories such as scalar-tensor gravity theory (27| 
and "brane world scenario" [28J predict large enhance- 
ment of the gravitational constant (in other words the 
reduction of 4 dim. Planck mass). There may have been 
an epoch in the early universe and may observe at an 
ultra high energy experiment where the gravitational ef- 
fects on hadrons are crucial. 



V. CONCLUSION 

In this Letter, we have studied the gravitational ef- 
fect to the dibaryons in the axially symmetric ES model. 
In particular, we have investigated gravity coupling con- 
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